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Abstract 



^ ' We propose a (theoretical ;-) model for quantum computation where the result can 

lO ■ be read out from the time average of the Hamiltonian dynamics of a 2-dimensional 

^^ I crystal on a cylinder. The Hamiltonian is a spatially local interaction among Wigner- 

V^ ' Seitz cells containing 6 qubits. The quantum circuit that is simulated is specified by 

CD ■ the initialization of program qubits. As in Margolus' Hamiltonian cellular automa- 

ton (implementing classical circuits), a propagating wave in a clock register controls 
asynchronously the application of the gates. However, in our approach all required ini- 
O ; tializations are basis states. After a while the synchronizing wave is essentially spread 

around the whole crystal. The circuit is designed such that the result is available with 
Ch I probability about 1/4 despite of the completely undefined computation step. This 

^ ' model reduces quantum computing to preparing basis states for some qubits, waiting, 

O^. and measuring in the computational basis. Even though it may be unlikely to find our 

specific Hamiltonian in real solids, it is possible that also more natural interactions 
allow ergodic quantum computing. 



1 Introduction 

The question which control operations are necessary to achieve universal quantum 
computing is essential for quantum computing research. The standard model of quan- 
tum computation requires (1) preparation of basis states, (2) implementation of single 
and two-qubit gates and (3) single-qubit measurements in the computational basis. 
Meanwhile there are many proposals that reduce or modify the set of necessary control 
operations (see e.g. ^ HI EI ,4, ^ ) . Common to all those models is that the program is 
encoded in a sequence of control operations. 

Here we consider a model which requires no control operations during the computa- 
tion since the computation is carried out by the autonomous time evolution of a fixed 



* email: {janzing, wocjan}@ira. uka.de 



Hamiltonian. The idea to consider theoretical models of computers which consist of a 
single Hamiltonian can already be found in [SIIZIIH]- However, these models are not 
explicitly designed for implementing quantum algorithms. We start from Margolus' 
approach since it has the attractive property that the Hamiltonian is a homogeneous 
spatially local interaction between cells of a 2-dimensional lattice and is therefore "rel- 
atively close" to interactions in crystals. Margolus' Hamiltonian implements the dy- 
namics of a classically universal cellular automaton (CA). In his two-dimensional model 
the front of a spin wave propagates in one direction over the surface and controls the 
updating of the cells. Even though there is no globally controlled clocking of the up- 
dates his local synchronization ensures that each cell is not updated until all relevant 
neighbors are already updated. In the Margolus scheme the computer is always in a 
superposition of many computation steps. At the beginning one has to prepare the 
wave front such that it mainly propagates in the forward direction. Such a state is 
not a computational basis state. We found it intriguing to use only basis states. Our 
goal was to reduce the required control operations to the absolute minimum: input of 
the initial state, the writing of the program and the readout of the classical output. 
The basis states we start with consist of components propagating forward and com- 
ponents propagating backward. Our circuit is designed such that even the backward 
computation leads to the correct result. When the time average of an appropriate ini- 
tial state subjected to the Hamiltonian dynamics is measured one obtains the correct 
result with high probability. The state tells us whether the result has to be rejected. 
Hence one may consider the procedure as a Las Vegas algorithm. Our Hamiltonian is 
a sum of operators which act on 10 qubits in contrast to the 2-dimensional Margolus 
cellular automaton which needs interactions between 8 qubits for universal classical 
computation. In 9^ one finds a 4-local Hamiltonian where the time average of a single 
qubit encodes the answer of a PSPACE hard problem. But the Hamiltonian has to be 
constructed for the specific PSPACE problem. The Hamiltonian is not homogeneous 
and is not appropriate for universal computation. 

The structure of the paper is as follows. In Section|21we choose a set of four 2-qubit 
gates which is universal for quantum computing. In Section |31 we construct a 4-qubit 
gate which includes all these 4 gates into one controlled gate. This makes the computer 
programmable. Then we describe how the synchronization scheme of Margolus is used: 
A wave front of a clock register propagating around the cylinder ensures that the 
programmable gates are applied in correct time order. This propagation is done by the 
evolution of an appropriate Hamiltonian. In Section |1] we describe the symmetry of 
the crystal by the crystallographic concept of Wigner-Seitz cells. In Section El we prove 
that the time average leads to the correct result. The readout of this result is explicitly 
described in Section |21 In Section [7] we briefly show that ergodic quantum computing 
can in principle solve all problems in polynomial space for all problems where usual 
quantum algorithms need only polynomial space. At first sight, this seems to be in 
contradiction to the fact that time steps of usual algorithms are translated to spatial 
propagation (as in P). 



2 Universal set of gates 

We recall |1U| that the following types of gates are sufficient for universal quantum 
computation. Let (C^)*^" be the state space of a quantum register. Then we consider 
the following two-qubit and single-qubit gates which are assumed to be available for 
every pair of qubits or every single qubit, respectively: 



1. The Hadamard gate on a single qubit: 

H : 



1 /I 1 



2. The controlled-phase gate 



+ |0)(0| ®i 



where |0), |1) are the canonical basis states of C and 1 is the identity. 

Note that an exact implementation of the SWAP gate is possible. Therefore, without 
losing universality, we allow the application of controlled phase gates only on adjacent 
qubits. 

We assume that gates acting on disjoint sets qubits can be applied at the same 
time step. We call such a time step a layer of the quantum circuit. The depth of the 
quantum circuit is the number of time steps. 

For reasons that shall be clear later we consider circuits U which have a special 
layer structure (see Fig. ^. Each time step consists of several gates with the following 
restrictions: 

• In even time-steps we allow only two-qubit gates acting on the qubit pairs (fc, k+1) 
with even k. 

• In odd time-steps we have only two-qubits gates on (A:, k + 1) with odd k. 
In this scheme we distinguish formally among four 2-qubit gates: 

t/oo:=l®l, Uqi:=1®H, Uw := H <S, 1 , Un := K{al/^) . (1) 

Using these gates, we construct a circuit U with the following properties: Let 
/ : {0, 1}" -^ {0, 1}*" be the function we would like to compute. The unitary U acts on 
the input, the output register, and some ancilla register and computes / in the sense 

U{\x) ®\y)®\Q... 0)) = |x) 8) |y © f{x)) |0 . . . 0) , 

where ® denotes the bitwise XOR. By construction, we have 

U^{\x) |y) |0 . . . 0)) = \x) © ly) ® |0 . . . 0) . 

Without loss of generality, we may assume that f{x) ^ for all inputs x by extending 
/ with an additional bit which is always 1. 

Note that there are quantum algorithms where /(x) is only computed probabilisti- 
cally. We will neglect this fact since it is irrelevant for the principles of our construction 
and would make the discussion unnecessarily technical. 
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Figure 1: Decomposition of arbitrary quantum circuits into layers of two-qubit gates [/qo, • • • , Un 
acting on adjacent pairs. 

3 Constructing the crystal Hamiltonian 

Usually a quantum circuit is considered as a sequence of gates. However, the usual way 
of drawing it (like in Fig. ^) suggests spatial propagation. Now we consider quantum 
circuits where quantum information is really spatially propagated and the time-axis is 
represented by the second dimension. 

Our circuit is wrapped around a cylinder. The cylinder is covered by c x /i squares 
("cells") of equal size. We have h (for "height") columns and c (for "circumference") 
rows. We need c > 2h for reasons which will be clear in Section |21 The columns 
correspond to the qubits of the original circuit and the rows to its time steps (see 



figH- 



Each cell (j, k) contains a data qubit. They form the data space 



In the j-th time step we apply all gates of layer j. A gate of the original circuit acting 
on the qubit pair {k, A; + 1) in level j translates to a gates acting on data qubits in cells 
(j) k), {j, k + 1), {j + 1, k), (j + 1, fe + 1). It applies the original two-qubit gate to the 
qubits in row j and propagates the information to row j + 1- Furthermore, the vertices 
between those 4 cells contain two program qubits which specify which one of the two- 
qubit gates in eq. (^ should be applied. Explicitly, there are two qubits between cell 
(j, k) and {j + 1, A: + 1) if both k and j are even or both are odd (see Fig. |2). For each 
vertex with program qubits we define the gate 



V:=W Yl ^'™ ® ^'- 

l,m€{0,l}x{0,l} 



(2) 



where Pjm '■= \jm){lm\ projects onto the state \lm) of the two-qubit program register 
at a certain vertex. W is the swap gate which exchanges the state of the qubit pairs 
(j, k) and (j + 1, k) and the pairs {j, A; + 1) and (j + 1, fc + 1). 





Figure 2: (Left) Cylindric crystal consisting of c x h cells. A pair of program qubits is located at 
the red points. The lines indicate the boundaries of a cell. (Right) The circuit wrapped around the 
cylinder. Every time when a two-qubit gate is applied the information of both qubits is propagated 
one row upwards. The output region consists only of trivial gates, i.e. the information is only 
propagated. 



Figure 3: Correct (left) and incorrect (right) walls. Putting a brick at position k,j corresponds to 
carrying out the gate in level j acting on the qubit pair (j, k) and (j, k + 1). 

This makes our system programmable and will be essential for achieving our goal to 
construct a universal Hamiltonian which can simulate all circuits. We will only write a 
program on some part of the cylinder because we need the other part as output region 
(see Sectional). As we have already stated, a computation would consist of applying 
all gates in row j in the j-th step. However, this requirement is unnecessarily strong. 
Actually, the only rule is that each gate in row j can only be applied if both gates 
in row j — 1 which contribute to its input have been applied. These synchronization 
rules can be visualized by building walls with bricks (see fig. El). The synchronization 
conditions mean intuitively that incorrect walls are not allowed. In order to make this 
analogy perfect we introduce dummy single qubit gates at the boundaries of odd rows. 

We would like to construct a Hamiltonian such that its autonomous time-evolution 
corresponds to a computation which respects these synchronization rules. Margolus 
|llj solved this problem by introducing clock qubits as follows. 

Each cell contains a clock qubit. Let Tic = (C^)*^ be the Hilbert space of all these 
clock qubits. Define the operator 



G 



j,k 



a'' O a'^ 
a (d a 



where the annihilation operators a act on the qubits (j, k) and (j, k+1) and the creation 
operators a^ act on the qubits (j + 1, k) and (j + 1, k+1). These operator Gj^k propagate 
two I's in the qubits (j, k) and (j, k + 1) one row upwards 
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where the left lower corner of the cell is at position (j, k). All other configurations are 
mapped onto the zero vector. Now we define the operator 



G 



j,k 



where j = 0, . . . , c 



1 and 

k = 



0,2, 
1,3, 



,h — 2 for j even 
,h — 3 for j odd 



(3) 



In contrast to Margolus we do not consider a cyclic system in both axis but only cyclic 
in one direction. This is because we think that a crystal with 2-dimensional torus 
symmetry seems less realistic. 
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Figure 4: Visualization of the basis states of Tic (the space of allowed clock states) as brick walls. 



At the boundary we define a family of operators which act on only two adjacent 
cells: For each odd j we set 



G 



i-i 



1 

V2 



G 



j,h-l 



1 

72 



t 



a 
a 



(4) 



where the annihilation operators act on the qubits (j, 0) and (j, h—1) and the creation 
operators act on the qubits {j + 1, 0) and {j + 1, /i — 1), respectively. These operators 
propagate a 1 in the qubits (j, 0) and {j, h — 1), respectively, one row upwards 
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where the left lower corner of the rectangle is at position (j, k) with j odd and k = 
0,h — 1. All other configurations are mapped onto the zero vector. Now we include 
these operators in the operator G. 

Now we define a G-invariant subspace Tie, interpreted as the space of correct syn- 
chronizations. Intuitively, it is spanned by the set of all basis states corresponding 
to correct walls. The position of the uppermost brick in each column is denoted by 
symbols 1 as in fig. [l] 

Lemma 1 (Synchronization space) 

Let Tic be the space spanned by those basis vectors \a), where a is a 0-1-matrix of size 
c X h satisfying the following conditions: 

(1) Each column contains a single 1, the remaining entries are all 0. 

(2) Let jk be the index of the symbol 1 in column k. Then for the indices of any two 
adjacent columns k and k + 1 we have 



(3) If k + jk is even then 
If k + jk is odd then 



\jk - jk+i\ < 1- 

jk > jk+i ■ 
jk < jk+i ■ 



Then Tic is G-invariant. 

Proof. Let a be any configuration satisfying the above conditions. 



* 








* 


1 





* 





1 









* 


1 


1 


* 








* 



Figure 5: Left: the application of the operators corresponding to the squares annihilate the state. 
Right: The application of the operator which corresponds to the upper square propagates both 
symbols 1. 

(1) Applying Gj^k to \a) does not lead to the zero vector iff the symbol 1 is at position 
j in the adjacent columns k and k + 1, i.e., jk = jk+i = j- If this is the case, 
then Gkj propagates both I's one position upward. Therefore, the configuration 
Gkj\a) still fulfills condition (1). 

(2) Assume first that a is a configuration with jk > jk+i for some k. Since a satisfies 
condition (2) we know that j^ = jk+i + 1. The only operators which act on qubit 
{k,jk) are Gj^^k and Gj^-i^k-i- The Gj^,^k operator vanishes when applied to \a) 
because the symbol 1 is at position jk — 1 in column k + 1 and not at position jk 
which would be required for a non-trivial action of Gj^k- The operator Gj^^i^k-i 
vanishes because the 1 is at position jk in column k and not at position jk — 1 as 
would be required for a non-trivial action of Gj^^i^k-i- This situation is shown 
on the left in fig. El 

The case jk < jk+i is proved analogously. 

(3) Assume that a is a configuration with jk = jk+i for some pair of adjacent columns 
k and k + 1. Let k + jk be even. In this case we have to show that every action 
which increases jk+i also increases jk- 

We first consider the case that both k and jk are even. The only operators that 
act on the qubit {k + l,jk+i) = {k + 1, Jfc) are Gkj^. and Gk+ij,.-i- The second 
operator vanishes because the symbol 1 is not at position position jk — 1 in the 
column k + 1 but at jk- The operator Gkj^ increases jk and jk+i as claimed. 
This situation is shown on the right in fig. [51 

Analogously, we can prove that this is also true if k and jk are both odd. The 
remaining case is that k + jk is odd. By using analogous arguments we can show 
that every action which increases jk also increases jk+i- 

D 



In analogy to Margolus' and Feynman's ideas we define the forward time operator F 

by 



F :-- 



j,k 



G,.k ® V, 



j,k 



where Vj^k is the gate V in eq. (O acting on the qubits {j,k),{j,k + l),(j + l,k), and 
(j + l,k + 1). For the operators Gj^^i and Gj^h~i (j odd) at the boundary we set 



Vj-i := 1 and Vj^h-i 



1. 



The Hamiltonian is defined as the sum of the forward time operator and backward 
time operator 

H ■.= F + F^ . 

In the sense of J12 l this is a 10-local interaction since each operator Gj^k ^-cts on 10 
qubits at once. Note that one may rewrite the interactions as /c-local terms with k < W 
by introducing qudits, i.e., particles with higher dimensional Hilbert spaces. Therefore 
the size 10 does not necessarily mean that this interaction is unphysical. 

To analyze the dynamical evolution we need the feature that F is a normal operator 
on the relevant subspace (analog to Margolus' results). However, since we do not work 
with cyclic boundary conditions the proof is a little bit more technical. As noted in jj^ 
the dynamics of the 1-dimensional cyclic Margolus Hamiltonian ^J| is the quasi-free 
time evolution of independent fermions. 

Even though we do not see if the clock dynamics of our Hamiltonian is also quasi- 
free, we can prove: 

Lemma 2 The restriction of F to the relevant space 

n := He (^ no (E) Tip 

is normal, i.e., FF^ = F^F. 

Proof. For an initial state \tlj) G TC the operator FF^t/j) is a sum of terms of the form 

FJ,kFl-J^P) . (5) 

F^F\ip) is a sum of terms of the form 

Flf,Fj,k\i^) ■ (6) 

For \k — k\ > 2 or |j — j| > 2 the operators F-~ and Fj^j^ act on disjoint qubits and 
thus commute. Then the products in eq. © and eq. @ are equal. 

li \k — k\ < 1 and \j — j\ < 1 then it is easily checked that the product G- rGj^k\a) 

is only non-zero for (j, k) = (j, k). 

Therefore, it is sufficient to show that 

(i,fc) U,k) 

in order to prove that F is normal. Since the operators Uj^^ are unitary it is sufficient 
to show that 

E^^>4fcl«) = E4fcG',,,|a) (8) 

ij,k) {3,k) 

for every allowed clock configuration a. Note that \a) is an eigenvector of the operators 
on both sides since each term which does not vanish is identical to a multiple of the 
vector I a). First we consider only the operators Gj^k which act on 4 clock qubits and 
not the special operators Gj-i and Gj^h-i ^it the boundaries. In the cyclic model of 
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Margolus, the right-hand term in © counts the possibihties to go forward and the 
left-hand term the possibihties to go backward. The fact that both numbers coincide 
prove normahty. in the non-cychc case the possibihties to add or remove half bricks 
have to be considered separately. 

Note that Gj^k\o) can only be non-zero if j = jk (with jk defined as in Lemma ^, 
i.e., there is the symbol 1 in position j in the fcth column. Then the term Gjj^^kW) does 
not vanish if and only if there is also a symbol 1 in position {jk, k + 1). 

To formalize these conditions we introduce the variable Ck '■= {jk + k) (mod 2) 
indicating whether jk + A: is even or odd. Due to the definition of the operators Gj^k 
the term Gj^^k\o) is only non-zero if Ck = 0. The position of the second symbol 1 
requires c^+i = 1. Since jk+i can differ from jk by at most 1 the conditions 

Cfc = and Ck+i = 1 

are also sufficient that Gj^,^k\o) is non-zero. 

Similarly, we have Gj^^^kGj ^|a) = \a) whenever 

Cfc = 1 and Ck+i = . 

Let nio and noi be the number of occurrences of the patterns 10 and 01 in the string 
(co, . . . , Ch-i), respectively. If nio = riQi then the leftmost and the rightmost symbols 
coincide. In both cases exactly one of the boundary terms 

does not vanish and yields the vector (l/2)|a). The same is true for the terms with 
the conjugated boundary operators. Hence both sides of eq. (jH| yield the same vector 
(n + l/2)|a). 

Note that tliq and tiqi can differ by at most one. This is the case if and only if 
the leftmost and the rightmost symbol are different. If nio = ^oi + 1 the leftmost and 
the rightmost symbols are 1 and 0, respectively. If ngi = nio + 1 they are and 1, 
respectively. In the first case only the combinations 



Gjo,-iGjo-i\a) , G^-^_^^;j_iGj^_i,/i-i|a) 



lead to non-zero terms and contribute to the right-hand side of eq. (jHI with (1/2) | a) 
each. The conjugated boundary operators lead both to vanishing terms. This fact 
compensates the difference of 1 = nio ~ ^oi in the contribution to the left-hand and 
the right-hand side of eq. 0. The second case (noi = nio + 1) is treated analogously. 

D 

The fact that F is normal helps to understand the dynamical evolution according 
to F + F^ . In jllj this fact makes it possible to find a conserved quantity interpreted 
as the computation speed. It is given by the operator V := {F — F^)/i. Then Feynman 
and Margolus start with initial states which have a positive expectation value of the 
computation speed. Their initial states are necessarily superpositions of basis states 
because the expectation value of V is zero for every basis state of the clock. Since \ip) 
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Figure 6: The black square is a Wigner-Seitz cell containing 6 qubits. The program bits are green, 
the clock qubits are red and the data qubits blue. The red grid indicates the original cells. 

is orthogonal to F\il)) and F^\il)) we have {ip\F — F^\ip) = with \ip) := |o) \(f)), where 
|i?i>) G Tip ® Tin and a is an allowed clock configuration. In our approach, all initial 
states are basis states. Despite these differences, normality of F will be essential in 
Section IS] for the "ergodic theory" of our Hamiltonian. 



4 Symmetry of the crystal 



The symmetry of a crystal can be described by a unit cell such that the whole lattice 
consists of shifted unit cells where the translations are integer multiples of the lattice 
vectors. A usual way to choose unit cells is given by the so-called Wigner-Seitz cell T^ . 
It is constructed as follows. Consider an arbitrary point Q in the lattice and consider 
the set of all points Q' which are equivalent to Q in the sense that the translation 
QQ' is a symmetry operation. Consider the perpendicular bisector of the side QQ'. It 
divides M^ into two half-planes containing Q and Q', respectively. Then the Wigner- 
Seitz cell (WS cell) is the intersection of all half planes containing Q. Here we choose 
the position of one pair of program qubits (i.e., a thick red point in fig. |2)) as Q. In the 
sequel we will refer to our original cells simply as cells in contrast to the WS cell. A 
WS cell is a square which has double area compared to the original cells and is rotated 
by 45°. It covers 4 adjacent cells such that it contains half of the area of each. This is 
depicted in fig. El 

We locate the data and clock spins of each cell such that the WS cell (which is 
centered around 2 program qubits) contains 2 clock qubits and 2 data qubits. Hence, 
the WS cell contains 6 qubits. Each WS cell interacts with those adjacent WS cell 
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which have an edge in common. Note, however, that the interaction among the WS 
cells are not pair-interactions between adjacent WS cells because it contains operators 
which act on 5 WS cells at once (note that the operators Fj^^ involve 4 cells). 

Note that the crystal is symmetric under reflections at columns, due to the sym- 

1 /2 

metry of the controlled A{az )-gate. The crystal, as we defined it, is not symmetric 
under reflections at rows. 



5 Mixing properties of the time evolution 

Our crystal Hamiltonian H = F + F^ is (on the relevant subspace) two times the 
real part of the normal operator F. Therefore, F and H have a common spectral 
decomposition. The following property of F is essential: 

Let \ijj) := I a) 10) € W be an initial state of clock, program, and data register, 
where a G {0, 1}'^'^ is an allowed clock configuration. Let g := {h + l)c/2 be the number 
of bricks (where half bricks are counted like full bricks) needed to cover the whole 
cylindric surface. Then we have 

F^lij) ±F''\iP) (9) 

for all j ^ k mod g. This is easily checked because each state F^\ijj) is a superposition 
of states where "the wall" is enlarged by j bricks. In order to get the same clock 
configuration one needs to add a multiple of g bricks. Note that the quantum circuit U 
(which is encoded in the program register) can be constructed in such a way that the 
orthogonality relation Q holds even for all j ^ k mod 2g. Consider, for instance, the 
case j = k + g. Project both states in Q onto the subspace of TC induced by a definite 
clock configuration a. On this subspace, the states of the data register differ by some 
unitary. This unitary U' is given by the concatenation of all those gates which have 
to be applied in order to go from the clock state a to a again by winding around the 
cylinder once. In other words, U' is obtained by splitting the circuit U and reversing 
the order of both parts as follows: Let Ui be given by the sequence of gates which are 
applied when the clock wave moves from its initial position to a. Analogously, U2 is 
given by all gates that are applied when the clock wave moves from a to the initial 
position. Then U' is given hy U' := U1U2 and U = U2U1. 

If at least one bit of the computed value f{x) is 1 the application of U leads 
always to orthogonal states in the data register^. Hence the orthogonality relation Q 
is already satisfied for j = g and k = 0. This corresponds to the trivial splitting 
Ui = U and C/2 = 1- Unfortunately, the bit flip which occurs on one of the output 
bits cannot be implemented by one gate since classical gates are not available in our 
setting. Therefore the other splittings may divide the flip operation into non-classical 
operations. In order to guarantee that also the other splittings lead to orthogonal data 
states we may construct U in such way that it flips two bits, one at the beginning and 
one at the end. Then either Ui or U2 contain one complete bit flip. 

The following lemma is important for analyzing the ergodic behavior since it shows 
that F is essentially a copy of the shift operator acting on mutually orthogonal spaces: 



^If the output part is correctly initialized. 
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Lemma 3 Let B he a normal operator on a finite- dimensional Hilhert space 7i and 
\^) GTC such that 



BH 



±B* 



j ^ k mod N 



for some A^ € N. Define 



H, := span,.eNo{i3'+''^l^)} 



,Af-l-- 



where I = 0, . . . ,N — 1 and Ti := ©;=o "^'^ Then only the following two cases can 
occur: 

1. All Til have the same dimension r. Then we may identify Ti with C ® C^ such 
that Til corresponds to |/) C and B (restricted to Ti) has the form 

B = S®A + S^®A^, 

where S is the cyclic shift on C^ and A is some normal matrix of size r x r. 

2. AllTii except forTio have the same dimensionr andTio has dimensionr+1. Then 
we may identify Ti with C © (C © C) such that Tio corresponds to C © (|0) © C) 
and Til to \l) fS> C for I = 1, . . . ,N — 1. Furthermore, this identification can be 
chosen such that the restriction of B to Ti. has the form 

B = Oe{S^A + S^0A^). 



Proof. Obviously B has the form 



/ 

^0 



An-1 \ 



>Aj^_iA]\f^i 



V An-2 / 

where each Ai maps from Tii to Ti^,-^ mod N' "^^^ diagonal entries of B^ B are 

The diagonal entries of BB' are 

AjV-l^Jv_l , AqAI ,... , ^^_2^jv_2 • 

Since B is normal we conclude 

ij+i) mod TV ^0+1) mod N - ^J^j 



Note that AjA^- and A^-Aj have the same rank, namely rank{Aj) = rank{A'-) 



(10) 



This 



AjTij for 



shows that all Aj have the same rank r. By definition we have Tij+i 
j = 0, . . . ,N — 2. Therefore, the dimension of Tij for j = 1, . . . , A^ — 1 is r. Only the 
dimension of Tio is not yet determined. Note that the dimension of Tio can not be 
smaller than the dimension of Tii (since the latter is the image of Tio)- 
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If B has only trivial kernel in Ti then Aq has also trivial kernel. Then the dimension 
of TCo is also r. This corresponds to the first case. 

The following arguments show that we can find a transformation which changes 
every Aj to the same matrix A. 

Let B = \B\U be the polar decomposition of B. Here U is unitary and \B\ := 
VBB'^. Note that U and \B\ commute since B is normal. Furthermore, \B\ has full 
rank and leaves each Tii invariant. 

We have UTii = ^-(i+i) mod N- Therefore, the power [/^ leaves each subspace 
Til invariant. Let X := vU^ in the sense that X is a function of U and not an 
arbitrary operator V with V^ = U^ . It also leaves each Tii invariant. Define A 
as the restriction of \B\X to TLq. We identify the subspaces Tii with ^((+1) mod 7V 
with each other via the unitary transformation X~^U. Note that X commutes with 
U. Therefore, this identification is consistent due to (X^^U) = X~ U = 1. By 
applying the transformation (X^^C/)^^ / times we transport the vector l'^) to the 
subspace TIq and obtain {XU~^y\Tp) . By applying A to this vector and transporting 
it back from TIq to Wj+i we obtain: 

{X-^U)^+^ A{XU-^)^\ij) = {X-^U)^+^\B\X{XU-^)^\iIj) 

= \B\XX~^U\ij) 

= \B\U\^l^) 

= B\^). 

This shows that our identification of subspaces allows to describe B by the action of 
the same operator A for each pair TCi and ^((+1) mod N- By choosing an arbitrary basis 
for Tio we may identify all spaces with C. This concludes the proof of the first case. 
If B has a non-trivial kernel in Ti. it is easy to see that its dimension is 1. This is 
due to the fact that the vectors B^^ \iIj) are in the image of B for all i > 1 and are 
orthogonal to its kernel. Then we may restrict B to the orthogonal complement of its 
kernel and obtain the first case. □ 

With the isomorphism of Lemma El we find statements about the time-average: 

Lemma 4 We adopt all notations of Lemma\^ For |i/j) G TCq define the time-average 
\t by 



^ Jt=o 
Let W be the probability measure on 0, . . . ,N — 1 defined by 

VF(/):=tr((|/)(/|®l)|V)(V|T). 

Let A = ^jttjQj be the spectral decomposition of A. Assume that {ip) lives in the 
subspace of eigenvalues of A with large modulus, i.e., \iIj) E X^j(1 ® Qj)'Hq where j 
runs over all indices with \aj\ > e. For 

16N 2N 
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with A := sin(((5/7)^) the total variation distance between W and the uniform distri- 
bution is at most 5, i.e., 

I 
Proof. We have |^) = |0) |k) with |k) G C^. Hence we have 

l^(/)=tr((|/)(/|0l)(|O)(O|0|K)(K|)T), 
where the time average is computed according to the Hamiltonian B-\-B'^ = S(^A-\- 

The projections 1 ® Qj commute clearly with \l){l\ ® 1 and with the Hamiltonian. 
Hence we can equivalently consider the time average of the mixture 

^(1 Qj)|0)(0| \ti){n\{l Qj) ■ 

3 

This is also true if we use 1-dimensional projections Qj instead of the original ones. 
On the image of each 1 ® Qj the time average problem reduces to the following 1- 
dimensional continuous quantum random walk according to the Hamiltonian 

Hj = QjS + TIjS' . 

where S is the cyclic shift on C^. Calculations on the explicit dynamics can be found 
in ^2l (for Uj = 1), we are only interested in time averages. We modify the techniques 
from J7] for studying discrete quantum walks to the continuous case. 

Let us now consider a fixed index j which is dropped in the sequel. Then we 
compute the probability distribution on 0, . . . , A^ — 1 induced by the time average of 
the state (|0)(0|)t according the Hamiltonian H := aS+aS^. Let a = re**^ be the polar 
decomposition of a. Then the eigenvalues of H are 2r cos(afc), where ak := <p + 2TTk/N 
for k = 0, . . . ,N — 1. Note that A^ is even. For simplicity we consider first the case 
r = 1/2 and derive an upper bound on the mixing time for this case. By rescaling the 
time we get a general bound. 

Now we consider the system with respect to the Fourier basis. Then we denote the 
eigenvectors of S with eigenvalue a;'^ := e^'"^^'^ by \k). The Fourier transform of the 
original basis states shall now be denoted by |e/). The initial state |eo) is an equally 
weighted superposition of all \k), i.e., the density matrix 

1 ^-^ 

We call eigenvalues and their eigenvectors good\i\sviiak\ > A with A = sin((5/7)^) and 
denote the number of good eigenvalues by A^'. The length of the interval of intervals for 
which I sin(a)| > A for a € [0, 27r) is 4arcsin A/(27r). Therefore, we have the following 
bound for large N: 

N-N' ^ 6arcsinA ^ 6(6/7)' ^ 5 , 

N - 2ti - 2tt -^7' ' ^ ' 
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Instead of the superposition of all eigenvectors we consider in the following an initial 
vector which is an equally weighted superposition of only "good" eigenvectors: 

|eo>-|/3>:=^ E 1^)- 
good k 

The trace norm distance between the modified density matrix and the true initial state 
is at most 

N -N' In -N' In -N' 

+ 2\ r^<3A 



N \l N - y N ' 

where the second term in the sum stems from dropping the bad eigenvalues and the 
first from rescaling the remaining part. Using eq. ()11|1 it is smaller than 

3(5 5 
— < - . 

7-2 

The time average of the modified state is 

PT := ^ r e-^^'\P){P\e^^Ut = ^^ T e^K'-"'=)*dt|A;)(A;'| . (12) 

^ -^0 ^^ ^ k,k' J^ 

The distance between two adjacent values aj- is 2tt/N . The derivative of the cosine is 
at least A or at most —A for good eigenvalues. Therefore, for a given k there is at most 
one k' such that | cos(afc) — cos(afc')l < Att/N, one in the interval where the cosine has 
negative derivative and one in the other interval with positive derivative. If we had 
three values a^, ak',ak" such that the distance between cos Uk and cos Ok' and between 
cos Ok' and cos ak" is less than An/N then we would have | cos ak — cos ak" \ < A2tt/N. 
Then we have at least two vales ak in the same interval which cannot be this close to 
each other due to the assumption on the derivative. 

Define projections Qp for every equivalence class p, i.e., Qp projects onto the span 
of all \k) with k G p. Note that these spaces are either 1-dimensional or 2-dimensional. 
We want to show that the probability distribution 

P{1) := {ei\pT\ei) (13) 

is almost the uniform distribution on the A^ points I = 0, . . . ,N — 1. We start by 
showing that the modified distribution 

R{1) := {ei\Y,QpPTQp\ei) (14) 

p 

is almost uniform. Explicitly, we have 



where 

fik,k') 



N'N 

k,k 



T lo e*^™' ""' '''°' ""'^^(^^ foi' ^ and k' equivalent 
otherwise 
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and the sum runs over all ordered pairs {k,k') of good indices. 

We measure the distance between the probability distributions P and R by the 
total variation distance 



p-m--=lY.\p{i)-m)\- 



2 I 

The Diaconis-Shahshahani bound ^H] estimates the total variation distance from R to 
the uniform distribution [/ by a sum over the Fourier coefficients of R: 

\\R-U\\<^^\R(m)\' 

Note that the first term of eq. H15|) has only a contribution to R{0). Hence we have 
only to consider the second term. We obtain 

7V-1 
1=0 k,k' 

where the last sum runs over k, k' such that k — k' = m mod N. There is at most one 
equivalent pair k, k' satisfying this condition. The reason is that one index k is in the 
region with negative derivative of the cosine and one index k' in the positive region. 
Let 1,1' be another equivalent pair where I is in the negative region. Since / and k are 
in the same region we may assume without loss of generality I = k + d with d < N/2. If 
/ — /' = k — k' mod A^ we must have I' = l + d mod A^. Then cos a/ < cos a^ — dA2Tr/N 
and cos ai' > cosa^' — dA2Tr/N. Hence / and /' cannot be equivalent. Therefore we 
find 

m^O {k,k') 

The last inequality is due the fact that there are at most A^' ordered equivalent pairs 
(of good eigenvalues). This proves ||-R — ?7|| < 1/A^ which is clearly smaller than (5/4 
for sufficiently large A^. 

Now we consider the total variation distance between P and R. Using the explicit 
representation ()12() of pT and the definitions of P{1) and R{1) in eq. H13() . ()14() we have 

\\p-m = ^Yl E K^H^)i K^'i^')i I £ e^(™^'^'=-'=°^"^')* dt\ , 

where the sum runs over all good inequivalent ordered pairs {k, k'). Note that we have 
|(e,|A;)| = l/^/N. Due to 

1;^ re'^'dt\<-^ 

T Jo '-r|x| 

we have 



N'T ^-^ cos ak — cos a^' 

k,k 
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For fixed value k we divide the inequivalent values cos(afc') in classes m = 1, . . . , [2A^/A] 

such that 

Am Aim + 1) 

— — < |cos(afc) -cos(afc')| < — . (16) 

The cosine function is on the interval [0, 2tt) two to one and its derivative has at least 
modulus A for the good eigenvalues. Therefore we have for a fixed k for every ttt, at 
most two k' such that 

Am Aim. + 1) 
cos Ofc' G [cos ak + -r— , cos ttk H — J , 

for 2N/A > m > —2N/A. Hence the inequality (|16jl is at most for 4 values fulfilled. 
Therefore have 

AN' ^^ 1 4 ,2N,^ ,2N, 

11^ -fill £ AFrrT,^,i^(^)M^y 

m=l 

In order to have this term less than 5/4 one has to wait the time 

327V , 2N 
A^^^X- 

Rescaling the dynamics with the modulus 2rj of the eigenvalues of A the time is 
increased by the factor l/min{2rj} = l/(2e). Then we obtain the time T as stated 
above. 

Putting everything together we obtain 

||VF-?7|| < - + \\P-R\\ + \\R-U\\ <6, 

where 6/2 stems from the restriction to good eigenvectors. □ 

For the initial state vector of our computation we have the problem that we do not 
know a priori whether a large component lies in the kernel of F. This is important since 
this component remains stationary under the evolution. For the part in the image we 
would like to know whether a large component lies in the subspace of small eigenvalues. 
This component would require a long mixing time. To address both problems we need 
the following two lemmas: 

Lemma 5 Let B be a normal operator on a Hilhert space and 1-0) an arbitrary unit 
vector. Let a be the angle between |^) and B\tjj) . Let P be the projection onto the 
image of B . Then 

tr(P|V')(V'|-P) > cos^a. 

Proof. The projection of |^) onto the image of B has at least the length of the 
projection of j-^) onto the span of -B|^) since the latter is a subspace of the image. 
Hence the projection onto the image has at least the length cos a. □ 

In order to estimate the mixing time we need the following: 

18 



Figure 7: Initial configuration where the component in the kernel of F can be estimated. The 
dashed brick indicates the only possibility to add a brick. 



Lemma 6 Let L := Hi?!";/;)!! be the length of -B|^) and < 5 < L. Let Ps be the 
projection onto all eigenspaces of F with eigenvalues of modulus at least 6. Then we 
have 

tr{Ps\ip){ip\Ps) > cos^ia - arcsin(VL)) , 

with a as in Lemma\^ 

Proof. Define the operator B^ := P^B. Due to \\B^ — L^\\ ^ ^ the tip of the vector 
B^\ip) is in an 5-sphere around the tip of B\il:). By elementary geometry, the angle 
between B^\'4)) and \il)) is at least a — arcsin((5/L). Since P^ is the projection onto the 
image of B^ we obtain the statement using Lemma [S] above. □ 

To use the lemmas above we could use the initial state a G {0, lY^ of the clock 
which is indicated by the wall in Fig[7| The only possibility to add a brick is at the 
rightmost position (cell l,h — 1). Hence F\a) = (l/v2)|a') where a' is the new wall 
with the additional half brick. To calculate F'\a') note that a' allows only two ways 
to remove a brick, namely that one just added (then a' is mapped to a again) and the 
upper most brick on the left (a' i-^- a"). This means that 

FtF|a) = i(|a)e|a")). 

Hence the angle between F'^Fla) and \a) is it/A. The length of -F^-F|a) is L = l/\/2. 
Using Lemma El we obtain 

tr{P\a){a\P) > cos2(7r/4) = 1/2, 

and by Lemma IHl we have 

tr{Ps\a){a\Ps) > cos2(7r/4 - arcsin(\/2<5) . 

Note that eigenvalues of F^F of modulus 6 correspond to eigenvalues of F with modulus 
y/S. If -Pf is the spectral projection of F for eigenvalues of at least modulus e we have 



tr{P,\a){a\P,) > cos^(7r/4 - arcsin(V2e^)) . 
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We will use the lemmas in this section to estimate the probability of success of the 
ergodic quantum algorithm. The idea is as follows. In Section |H1 we will argue that 
the correct result can be found for all states -F-'IV') (where \ip) is the initial state) for 
all j satisfying a certain condition. This condition ensures that the circuit U has been 
applied an odd number of times on the data qubits. To formalize this we introduce 
spaces Til ss in lemma IHl which are spanned by the vectors F^-' |^). Let Qi be the 
projection onto 7ii. Then for the probability distribution induced by the time average 

W{1) = ^ j^ tr{Qle^-^'\^l^){^l^y^')dt 

we have a lower bound on each W{1): 

Lemma 7 There is an initial state of the clock configuration such that the probability 
W{1) to find the state in TCi after one has waited the time T as in lemma^is at least 

W{1) > (^ - 2(5) cos2(^/4 - arcsin(\/2e2)) . 

The proof follows immediately from the lemmas of this Section: We choose the 
initial clock configuration of fig. [3 Above we have argued that the probability for 
finding a state in the eigenspace with eigenvalues of modulus at least e is given by 
the cos^-expression on the right. Given a state in this subspace we have uniform 
distribution up to a variation distance 5. This yields the factor 1/A^ — 26. 

In the next section we show for which part of the spaces TCi we have certainly a 
correct result and how this promise is used in the readout procedure. 

6 Initialization and Readout 

It is clear that the program qubits have to be initialized according to the simulated 
quantum circuit. Furthermore we have to initialize all clock qubits. On the data 
register we have only to initialize those qubits which are located in the cells where the 
initial clock wave is located. 

The readout of the computation result is done as follows. Here we assume that the 
initial state of the clock register is \a) where all symbols 1 are in row (In Sectional 
we have also considered another initial configuration which makes is easier to decide 
which component of the initial state is in the kernel of the Hamiltonian. However, the 
analysis of this Section is technically more complicated and the output region would 
have to be enlarged for this initial configuration). We define an output region which 
consists of all cells with column index between 1 and m where m is any natural number 
greater than 2h. 

We choose an arbitrary row in the output region and measure as many clock qubits 
of this row as are necessary to find the wave front. If we have found a clock qubit 
in state 1 1) in position j, k the wave front in row j + 1 and j — 1 has to be in one 
of the columns k — 1, k, or k + 1. By this procedure we can localize the whole wave 
front. If it is completely localized in the output region we know that the state of the 
corresponding logical qubits is either of the states U{\x) (gi |0 . . . 0)) or \x) (d \0 . . . 0) . 
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Then we can readout the result. We may define / in such a way that we can decide 
whether the result is correct or not. In the following we will give a lower bound on the 
success probability of the whole readout procedure. First we estimate the probability 
for finding the wave front in the output region. 

The wave front starts in row 0. States in Tii are in general superpositions of different 
wave fronts. Note that every such wave front consists of I mod g bricks. By elementary 
geometric arguments one can check the following statements: First we consider the case 
that I is in the interval 0, . . . ,g — 1. A wave front which consists of more than h'^/2 
bricks has completely passed row 0. Similarly, all row indices of the symbols 1 can be 
guaranteed to be at most m if / is at most h{m — h)/2. Therefore, we have at least 
h{m — h)/2 — h?'/2 = h{m — 2h)/2 spaces TCi which are completely in the output region. 
We obtain the same number of spaces Tii for I = g, . . . ,2g — 1. By these arguments we 
can easily derive the following lower bound. For each space Tii we can guarantee at 
least the probability 1/A^ — 25. This yields the following bound. 

Lemma 8 (Probability to find the wave front) 

The probability for finding the wavefront completely in the output region is at least 

him-2h){^-25)s,, 

where s^ is the size of the component of the initial state in the eigenspace of F with 
eigenvalues of modulus at least e. 

However, only the values above in the second interval ensure correct output. Since 
the function / is without loss of generality 1 on at least one bit, we can distinguish 
whether the result has to be rejected and the experiment has to be repeated. The 
probability that the first experiment succeeds can hence be estimated by dividing the 
lower bound of Lemma |H1 by two: 

Theorem 1 (One Shot Success probability) 

The probability for finding the output region and furthermore obtaining the correct 
computation result is at least 

2 ^N ' 

with Se as in Lemma\^ 

Note that we have chosen the initial state of Fig. [7| because we were able to prove 
a lower bound on the length of the component in the image of F which yields a good 
value for s^ ■ Actually, the disadvantage of this initial state is that the propagation is 
slow since the wall can grow only at one point. The more natural initial configuration 
given by a flat wall allows propagation in every second column. For this wall, we do 
not have a good estimation for the component in the image of F which is as simple as 
for the wall in Fig. Q Nevertheless, we believe that it is better to start with the flat 
wave. If the size of the output region dominates the size of the circuit (i.e. c^ m and 
c ^ h) the quotient h{m — 2h)/N tends to 1 since N = c{h + 1). With small e the 
factor in Theorem ^ is almost 1/2. Hence the success probability tends to 1/4. 
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7 Solving PSPACE problems in crystals of poly- 
nomial size 

It seems to be a general property of our construction that the size of the crystal 
necessarily grows linearly with the running time (i.e., the depth) of the encoded circuit. 
From the complexity theoretic point of view, this would have important consequences. 
Note that the complexity class PSPACE contains all problems which can be solved using 
polynomial space resources [W| . The running time of an algorithm solving a problem in 
PSPACE may be exponential. This seems to imply that the ergodic quantum computer 
would need exponential space in contrast to usual models of computation (e.g. Turing 
machines and Boolean circuits). Now we want to show briefly that even the ergodic 
quantum computer can solve all problems in PSPACE in polynomial space. 

The key idea is that even if an algorithm has exponential running time, it has nec- 
essarily (by definition) a polynomial description of the required sequence of operations. 
Therefore it is always possible to construct a circuit U of polynomial depth such that 
the repeated application of U solves the PSPACE problem. 

In f2D; we have shown that for every problem in PSPACE there is a two-gate 
quantum circuit U of polynomial size which computes a function / : {0, 1}" -^ {0, l}™ 
in the following sense: 

1. There is a (possibly exponentially large) natural number r such that 

C/'^(|x) |y) |0 . . . , 0)) = Ix) |y © f{x)) © |0 . . . 0) , 

where x £ {0, 1}" is the input string and y is an arbitrary string in the output 
register^ . 

2. The change of the state of the output register given by 

y^ y® f{x) 

occurs for a certain power s of U , i.e., for all U^ with < j < s the output state 
is still y and for all U^ with s < j < r — 1 it is already y © f{x). 

Furthermore, r and s are known by construction of U . This is possible since there is 
always an upper bound on the running time of an algorithm derived from the restricted 
space resources. By introducing idle cycles (counting steps) one can guarantee that 
this bound is exactly attained. Note that it does not make sense to require that the 
change of the output state occurs during the rth application of U . Otherwise / could 
be computed by a single application of C/^^. This is shown by the following argument: 

Assume 

U''{\x) © |0 ... 0) © \y)) = |x) © |0 . . . 0) © |y © f{x)) 

and 

C/'^-i(|x)©|0...0)©|y)) = |V)®|2/), 



^The construction in (20^ is restricted to binary functions. However, the generalization to several output 
qubits is straightforward. 
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where I?/') is an appropriate state of ancilla+input register. Then we have 
U-\\x)^\0...0)^\y®f{x))) = U-^U'''{\x)^\0...0)^\y)) 

This mean that one apphcation of U^^ maps y © f{x) onto y, i.e., is mapped onto 

The construction of 20 follows the usual philosophy of reversible computation |21j : 
The actual computation is done during the first r/2 cycles. Then the result is copied 
to the output register with Controlled-Not gates. The only goal of the last r/2 cycles 
is to undo the computation and restore the initial state. 

The ergodic theory in Section \^ applies directly to PSPACE problems after sub- 
stituting N = 2g to N := 2rg. Furthermore one has to guarantee the orthogonality 
condition Q for all j ^ k mod 2rg. The bit flips which have been explained at the 
beginning of Section El have to be substituted by incrementing counter registers. 

The readout is done exactly as in Section |H1 Given that we have localized the clock 
wave front in the output region we have the correct result with probability 1/2. As 
in Section El we can choose / in such a way that it indicates whether the result has 
to be rejected. Hence the probability of success is not reduced by the fact that the 
computation requires more cycles of U. 

8 Conclusions 

We have proposed a model of quantum computing which does not require any control 
operations during the computation process. The only required operations are the 
initialization of basis states and the readout in the computational basis. 

The relevance of this model is two-fold: first it shows that, in principle, quantum 
computation can be realized with a small amount of quantum control. Even though our 
interaction is rather artificially constructed, it is a priori not clear that it is unphysical: 
It consists of finite range interactions among cells of a crystal which contain some finite 
dimensional quantum systems. This shows that relatively simple local interactions in 
homogeneous solid states may have universal power for quantum computing without 
external control. We admit that it seems difficult to decide whether the interactions 
in real matter have such properties. However, this may be an interesting question for 
future research. 

The second aim of this paper concerns the thermodynamics of computation. As in 
IZIJH] the computation is performed in an energetically closed physical system with 
the additional feature that only the preparation of basis states is required. 

It would be desirable to find more simple Hamiltonians which are universal for er- 
godic quantum computing. A basis to find them could be given by simple 1-dimensional 
universal quantum cellular automata. 
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